Nonequilibrium Magnetization Dynamics of Nickel 
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Ultrafast magnetization dynamics of nickel has been studied for different degrees of electronic 
excitation, using pump-probe second-harmonic generation with 150fs/800 nm laser pulses of various 
fluences. Information about the electronic and magnetic response to laser irradiation is obtained 
from sums and differences of the SHG intensity for opposite magnetization directions. The classical 
M(T)-curve can be reproduced for delay times larger than the electron thermalization time of about 
280 fs, even when electrons and lattice have not reached thermal equilibrium. Further we show that 
the transient magnetization reaches its minimum as 50 fs before electron thermalization is completed. 

PACS numbers: 42.65.Ky , 75.40. Gb , 78.47.+p 



Ultrafast spin dynamics in ferromagnets is of great in- 
terest from both theoretical and experimental points of 
view. In particular, the short-time dynamics of mag- 
netism in transition metals, with many excited electrons 
^ • not at equilibrium with the lattice, is a new area of 
OO \ physics. Such studies are important for developing a 
theory of transient magnetization behavior in the sub- 
picosecond range. It seems that the only experimental 
data which can guide theoretical analysis are the ones 
t^- ' reported by Beaurepaire et al. jl| on time-resolved de- 
magnetization of Ni induced by femtosecond laser pulses 
of 620 nm at one specific fluence. The authors utilized the 
magneto-optical Kerr effect to detect hysteresis loops for 
different time delays between pump and probe pulses, 
i-^j ■ By comparing the time-dependent remanence with the 
equilibrium temperature dependence of magnetization, 
O ■ M(T), they derived the time evolution of the spin tem- 
perature within the framework of the phenomenological 
three-temperature model 0. Clearly, it is of great im- 
"kp | portance to confirm whether or not M(T) can be used 
to describe the transient magnetic response to electron 
excitations in itinerant ferromagnets and whether there 

_ _ i 

is a time delay between electron thermalization and mag- 
netization changes. 

In this Letter we present time-resolved data on 
the transient magnetization measured by pump-probe 
second- harmonic generation (SHG). The great advantage 
of this technique is that it allows to simultaneously fol- 
low electron-temperature relaxation and transient mag- 
netization, without further need for additional calibra- 
tion measurements. This is a consequence of the even 
and odd contributions to the nonlinear susceptibility || . 
The measurements were carried out for a large variety of 
pump fluences leading to different initial electron tem- 
peratures. After equilibration of the electron bath, we 
find the transient magnetization to be governed by the 
electron temperature T e via the classical M(T)-curve 
j|. However, we observe a strong deviation of the data 
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from the magnetization curve in the short-time range 
t < 0.3 ps, where the electron system is in a nonequi- 
librium state. 

When discussing nonequilibrium spin dynamics in fer- 
romagnets, it is worth noting that relaxation times pre- 
sented here, where nonequilibrium electronic states are 
excited, are not compatible with the spin-lattice relax- 
ation times, t s i. The latter denotes the time it takes for 
the spin system to adjust to a new temperature while 
equilibrium between electron and lattice temperature is 
maintained. Hence, the values of t s i = 100 ± 80 ps mea- 
sured by Vaterlaus et al. for ferromagnetic Gd |2||| and 
the limits given for Fe MM and Ni JjJ must be considered 
in a different context. 

Experiments were performed at room temperature on 
polished polycrystallinc nickel (purity 99.90 %) in air. 
The samples were magnetized to saturation with the 
magnetization direction oriented perpendicular to the 
plane of incidence. P-polarizcd 150fs/800nm pump and 
probe pulses with an intensity ratio of 3:1 for the highest 
pump fluence were used to excite the electrons and probe 
electron temperature and magnetization by detecting the 
reflected unpolarized SHG. Angles of incidence were 20° 
for the pump and 45° for the probe beams. We do not 
expect any interfering effects from the nickel oxide layer 
— which is antiferromagnetic — since we only measured 
relative changes in SHG yield for opposite magnetization 
directions j8| as a function of delay time between pump 
and probe pulses. 

Typical relaxation curves are shown in Fig.|l] for seven 
selected pump intensities, although more were measured 
and included in Fig. 3. Displayed are sums and differ- 
ences of SHG signals for opposite magnetization direc- 
tions / ± (t), normalized to the SHG yield I^ 1 for the probe 
pulse alone: 

A/ ± = (/ ± (t)-/ ± )// ± , (1) 

where 7 ± = 7(2w, +A?) ± J(2w, —M). The second-order 
polarization in the presence of magnetization can be writ- 
ten in the form || 

P l (2u;) = Y / [x^ n ( ± M) ± X^ d (±M)]E J E k . (2) 

The even tensor elements are almost unaffected by the 
magnetization, whereas the odd ones scale linearly with 
M Hence, we can substitute: 

evenfi\/f\ even odd/i\T\ odd fo\ 

Xijk ( M ) = Xo,ijk . Xijk ( M ) = Xo,ijk ■ \ M \ ■ I 3 ) 
From Eq. ^| we obtain the expressions 

/+ = 2./ 2 H. [\A X ^ en \ 2 + \B X ° dd M\ 2 ] , (4) 
r = 4 • J 2 H • \A X r n Bxt d M\ • cos^ , (5) 

where the individual tensor elements are combined to 
an effective even and odd contribution denoted by Xo ven 
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and Xo ■ A and B represent the corresponding effective 
Fresnel coefficients for the fundamental and frequency- 
doubled light 0. The resulting phase of A, B, xg" e ™, 
and Xo dd depends on M and is given by <j>. Measurements 
of transient linear reflectivities (not shown here) verify 
that the magnitude of the Fresnel factors A and B is 
nearly unaffected by T e . Also, calculations reveal that the 
magnetization-independent factors Xo" en an d Xo dd var Y 
only a few percent with T e for T e < Tq (Tq =Curie 
temperature) . This can be understood by the large width 
(« 0.3 eV) of vacant d-states in the minority band |ll| in 
comparison with the maximum electron temperature of 
550 K reached with the highest pump flucnce. Thus, we 
conclude that the time dependence of / is contained in 
M(t). Assuming M[t) = M(T e (f)), the classical M (T)- 
curve B can be approximated by M(T e ) = M(T )[1 - 
const. (T e — To)] 1 / 2 for the electron-temperature range 
covered in our experiment. Here, To denotes the mean 
electron temperature at negative delay. Inserting this 
into Eqs. (||), (||), and (j|) we obtain 

AI+(i) = const.- [T -T e (t)] , (6) 

arm = M{T f^ - 1 ■ P) 

Prerequisite for the validity of these equations is a ther- 
malized electron distribution. Thus for t > 0.3 ps the 
curves in Fig. la represent the time evolution of electron 
temperature, while those in Fig. lb describe the transient 
magnetization, provided <fi is constant. That A/ + is in- 
deed proportional to temperature can be verified by plot- 
ting the averaged values for delay times > 3 ps in Fig. la 
against flucnce. The result is shown in Fig. 2. In ther- 
mal equilibrium with the lattice, T e scales linearly with 
fluence. The slope m = 1 in Fig. 2 proves this proportion- 
ality in Eq.|^ and hence the validity of our approximation 
for M(T e ) for longer delay times. 

Eqs. |^ and [t] suggest a graph of (A/ _ + l) versus — A/+ 
which can be compared with the classical M(T)-curve 
for nickel Q. Such plot is shown in Fig. 3. In the upper 
frame [Fig. 3a], only data for delay times > 0.3 ps have 
been analyzed. Obviously, this plot reproduces within 
±5% the equilibrium magnetization curve, indicated by 
the solid line. The good agreement justifies the above 
made assumptions and demonstrates that (j> does not de- 
pend significantly on T e (t). Furthermore, it provides an 
intrinsic electron-temperature calibration for excitation 
with various pump intensities. It must be emphasized 
that this agreement includes the range between 0.3 ps 
and 3 ps where T e is not in equilibrium with the lattice 
temperature [cf. Fig. 1] . We also conclude from Fig. 3a 
that the electron temperature governs the magnetization. 
From the delay time where the points fall on the magne- 
tization curve we can unambiguously conclude that the 
electrons are fully thermalized after w 280 fs, when the 
minima in Fig. la are reached. 



3 



The nonequilibrium situation is demonstrated by in- 
corporating the data for t < 0.3 ps into the same plot, 
as illustrated in Fig. 3b for three fluences. The strong 
deviation of the data from the magnetization curve can 
be attributed to nonequilibrium electrons, and, in view 
of Eq. (||), to a rapid phase change. It vanishes when the 
thermalization of the electron gas is completed, and then 
the data again follow the magnetization curve. 

The different behavior of nonequilibrium electrons be- 
comes more evident in the real-time plot of A/ + and A/~ 
in Fig. 4, which is a magnified version of the short-time 
range in Fig. 1. For clarity, only data for three fluences 
are displayed, but the same delay between the minima of 
A/ + and A/~ is observed for all fluences. The minima 
of A/+ at (280 ± 30) fs (dotted line in Fig. 4a) mark the 
time at which the thermalization of the electrons is com- 
pleted and agrees with the value reported by Beaurepaire 
et al. . The breakdown of A/~ in Fig. 4b follows the 
convolution of pump and probe pulses. For all fluences 
the minimum is reached about 50 fs before the electron 
temperature is established. To understand the faster re- 
sponse of AI~ oc M compared to A/ + we notice that 
for t < 0.3 ps Eq. (||) is no longer valid and AI + oc |M| 2 . 
This means that A/ + reflects the momentary popula- 
tion difference between minority and majority spins and, 
therefore, monitors the development of electron temper- 
ature out of a highly nonequilibrium state ]I2] ]. A/~ is 
in addition sensitive to the direction of M and samples 
preferentially the minority spins of the vacant d-states. 
Note, that the spin-orbit coupling is restored faster for 
minority than for majority spins Jl3[ ] , since the width of 
the vacant d-states (« 0.3 eV) is large compared to the 
Fermi distribution for T e < Tc- The shift in Fig. 4 and 
the fact the magnetization curve is reproduced after ther- 
malization of the electrons prove that for t > 280 fs the 
electron temperature governs the magnetization. This is 
in striking contrast to the results reported by Beaurepaire 
et al. who observed a 2 ps delayed magnetic response 
with respect to T e . Further investigations are required to 
elucidate the cause of these deviating observations, since 
too many parameters differ in the experiments reported 
here and in [Q. Note, however, that the response time 
of itinerant magnetism is expected to be of the order of 
Tq , corresponding to about 80 fs for Ni. 

Let us return to Fig. 1 and attempt to understand 
the large decrease of AI^ in terms of the band struc- 
ture. The fundamental wavelength of 800 nm (1.55 eV) 
matches well the difference between the two mains peaks 
of the Ni minority band (jnj, as sketched in Fig. 5. Due 
to the exchange splitting, the upper peak is unoccupied 
for minority spins and can resonantly enhance the SHG 
in the intermediate state. For majority spins there is no 
such enhancement. The filling of the vacant states in the 
minority d-band causing the breakdown of M reduces 
this resonance enhancement which leads to the fiuence- 
dependent drop of signals in Fig. 1. The effect is further 
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amplified by the large matrix elements for d-d transitions 
JTij . The size of reduction of A/ + is determined by the 
ratio of odd to even parts of the nonlinear susceptibility 
and the corresponding Fresnel factors. Using the temper- 
ature calibration via Figs. 2 and 3, one can derive from 
Eqs.(3) and (4) the ratio \Bx Q dd \/\Ax%' en \ = 0.7 ±0.1. 

Summarizing, we have shown that pump-probe SHG 
is a powerful method for investigating ultrafast spin dy- 
namics in ferromagnets, since it contains information 
about both, electron temperature and magnetization. 
We applied different pump fluences to study the elec- 
tron and spin dynamics for various degrees of excitation. 
Independently of pump fluence, the electron thermaliza- 
tion is completed after (280±30) fs and the magnetization 
curve @] is reproduced from this time on, even when the 
electrons are not in equilibrium with the lattice. For the 
first time we have observed that the magnetic response 
is faster than the electron thermalization, which may be 
attributed to different lifetimes of excited majority and 
minority electrons. Further investigations of this topic, 
in particular frequency-dependent measurements which 
may uncover band-structure effects, are desirable. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft, Sfb 290. 
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FIG. 1. Time dependence of normalized sums (a) and 
differences (b) of SHG yields for opposite magnetization di- 
rections, as denned in Eq. (1). The curves were recorded with 
different relative ffuences, calibrated by 1.00 ~ 6 mj/cm 2 . 
Constant levels for t > 3ps reflect equilibrium between elec- 
tron and lattice temperatures, T e and T(. 
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FIG. 2. Normalized sums of SHG yields for opposite mag- 
netization directions, averaged over all values for t > 3 ps, as 
a function of relative fluence. The linear slope indicates that 
AJ + is proportional to electron temperature [see Eq. (6)]. 
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FIG. 3. (a) Comparison of all measured data for t > 0.3 ps 
with the equilibrium magnetization curve of Ref. [4], repre- 
sented by the solid line, (b) Data for three selected fluences 
covering the total time range. Deviations from the magneti- 
zation curve point to a nonequilibrium state of the electron 
and spin systems for t < 0.3 ps. 
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FIG. 4. Magnified version of Fig. 1 for three selected flu- 
ences. The minima at (280 ± 30) fs in (a), marked by the 
dotted line, signal the onset of a well defined electron tem- 
perature. The minima in (b) occur 50 fs faster (dashed line). 
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FIG. 5. Density of states for minority and majority spins 
according to Callaway and Wang [11]. While single-photon 
excitation proceeds with comparable rates from both parts, 
SHG is resonantly enhanced by the unoccupied peak in the 
minority band. 
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